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Abstract. In the 1990's exponential-type error bounds appeared in the 
theory of radial basis functions. For multiquadric interpolation it is 0{Xd) as 
d —>■ 0, where A is a constant satisfying < A < 1. For Gaussian interpola- 

tion it is 0{C'd)'d as d — > where C and c' are constants. In both cases the 
parameter d, called fill distance, measures the spacing of the points where 
interpolation occurs. This kind of error bounds is very powerful. However 
it only measures the difference between the approximant and approximand. 
Mathematicians and engineers often need to know the matching of the deriva- 
tives when dealing with partial differential equations. In this paper we extend 
this kind of error bounds to a form which measures the difference between 
the derivatives of the approximant and approximand. 
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1 Introduction 



In this paper h always denotes a continuous function on i?" which is condi- 
tionally positive definite of order m. The interpolation occurs at scattered 
data points {xj, yj), j = 1, . . . , N, where X = {xi, . . . ,xn} is a subset of i?" 
and the y'jS are real or complex numbers. We are going to use as interpolant 
function of the form 

TV 

s{x) = p{x) +'^Cjh{x — Xj), (1) 
i=i 

where p{x) is a polynomial in Pm-i and the c'jS are real or complex numbers 
satisfying 

N 

J2cjq{xj) = 

i=i 

for all polynomials q in Pm-i and 

N 

Pi^i) + Cjh{xi - Xj) =yi, i = l,...,N. (2) 
i=i 

Here Pm-i denotes the class of polynomials of i?" of degree less than or equal 
to m — 1. 

A sufficient condition for a unique solution of (1) is that h be strictly 
conditionally positive definite when X is a determining set for Pm-i- Recall 
that X is said to be a determining set for Pm-i if P is in Pm-i and p vanishes 
on X imply that p is identically zero. 

What's noteworthy is that in this paper both the interpolating and inter- 
polated functions belong to a function space which was also used by Madych 
and Nelson in [7] where they first introduced the exponential-type error 
bounds. We do not shrink or widen the function space. If h is continu- 
ous and conditionally positive definite of order m, then it induces a native 
space denoted by Ch,m which is a semi-Hilbert space for m > and a Hilbert 
space for m = 0. The definition and relevant properties of the native space 
can be found in [5] and [6J. It's not very easy to introduce this space thor- 
oughly in limited pages. However Luh made a lucid characterization for this 
space in [1] and [2]. For a better understanding of this space, we suggest that 
the reader read [T] and [2] first. 
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2 Multiquadric Interpolation 



Although Madych and Nelson made an incontrovertible contribution in the 
construction of the exponential-type error bounds, there are some drawbacks 
in their error bounds. For example, the calculation of the crucial constants 
in their error bounds is a big problem. Moreover, strictly speaking, their 
error bounds essentially apply to only. In order to make them suitable for 
i?", n > 1, some modification is necessary. Recently Luh made a thorough 
treatment for these problems in [3] and Of course it doesn't lower the 
author's respect for them. However, we choose to cite Luh's error bounds in 
[3] and [1] since they are more complete and it will make this paper more 
useful. 

For any set Vl C and X = {xi, . . . , , Xjy} C R^^ let 
d{Q,X) = supy^ninf^^x\y - x\ 

be the fill distance which will be denoted by d for simplicity in the following 
discussion. Now we need a lemma which is just Corollary2.5 of [3j with some 
notations simplified. 

Lemma 2.1 Suppose h is defined by 

h{x) := r(-|)(c2 + IxHf , (3 e i?\2iV>o, c > 0, 

where \x\ is the Euclidean norm of x in R^ , T is the classical gamma function 
and (3, c are real constants. Let m = [|] for /5 > and m = for (3 < be 
its order of conditional positive definiteness. Then, given a positive number 
bo, there are positive constants do and A, < A < 1, which depend on bo, c 
and 13 for which the following is true: If f & Ch,m o-nd s is the interpolant 
defined in (1) that interpolates f on a subset X of R^ , then 

\f{x)-s{x)\<C-{xt^-\\f\U (3) 

holds for < d < do and all x in a set Q C R^, where Q is any set in which 
for any point x, there exists a set E such that x ^ E O fl where E is a cube 
or the rotation of a cube of side bo satisfying the property that every subcube 
of E of side 2d contains a point of X. Here, \\f\\h denotes the h-norm of f 
in Ch,m, and C is a constant depending on n, c and p. 
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Remark. For the calculation of C and A in the preceding lemma we refer 
the reader to [5]. As for Ch^m and \\f\\h, we suggest [I] , [2] , [S] , [S] and [7] for 
detailed introduction. However, if the reader needs a cursory understanding 
only, maybe [7J is enough. 

Theorem 2.2 Let h, X and s be as in Lemma2.1. For any bo > 0, let do, A 
and Q be as in Lemma2.1. Let d = d{fl, X) be the fill distance, and m be the 
order of conditional positive definiteness of h. Suppose I > max{l,m} is an 
integer. Then 

(a) for every f G Ch,m and \a\ = I, the distribution D^f belongs 
to C{R^) and there is a constant such that for all f in Ch^m, 

\\Dy-D^sU<cJfU; (4) 

(b) for any y E Q such that the closed ball B(y,S) C n for some 
6 > 0, and any f e Ch^rn, 'if ^ < W\ < I, then 

\D^f{y) - D-s{y)\ < q^iMo"^ ■ Jlf (5) 

where C\a\ is a constant independent of f, Mq = C||/||/iA^ as 
m (3), and Ml = max{Mi, Mq/I^"') with Mi := C'\\f\\h where 

Proof. To prove (a), note that s G Ch^m jSj and hence / — s G Ch^m- Then (a) 
follows immediately from Proposition4.1. of [5] and the discussion preceding 
it. What's difficult is to show (b). 

For y e with 'B{y, S) C n for some 5 > 0, t e [-5, 5] and u e R"' with 
|m| = 1, let ipit) = f{y + tu) — s(y + tu). Then 

^(^){t) = k\ Y: ^D^{f-s){y + tu) 

\a\=k 

ior < k <l. By Proposition4. 1 of (6|, H^^'^IU < C'\\f - s\\h with C = 
^! Z]|a|=/ Ca/a!. The discussion of p. 223 of [6] tells us that ||/ — s\\h < ||/|U- 
Thus ||V'*-'''||oo < From (3) we also have a bound for \ip{t)\ on the 

interval [—6, 6]. By the result of Gorny in p. 139 of [8J we get for < < /, 

|^('=)(0)| < 16(2e)'^Mo"^ -Mj 
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where Mq = C\\f\\h\^ and M = max{Mi, MqI\5-^) with Mi = C'\\f\\h where 
C — l\ Y,\a\=i Ca/oi\. This imphes 



\a\=k 



< il6(2e)'=Mo -M} 



k fe 

i 



Since m is any unit vector in i?", we have 



SUP\u\=l 



7/" 

E -T^"(/-^)(2/) 



|a|=fe 



a! 



for every y & Q with B{y,5) C Q for some 5 > 0. Now, let T4 be the space 
of vectors v — {va)\a\=k- Then 



2^ 7T 



\a\=k 



is a norm for T4. Any two norms on a finite-dimensional space are equivalent. 
Therefore this norm is equivalent to the Euclidean norm. This gives that for 

each |a| = A;, 



|^"(/-s)(|/)|<^16(2e)^Mo-^M, 
for some constant pk > 0. Let = ^16(2e)'^. Then (b) follows. 



, 1_£ £ 



Remark: Note that when d is small, Mi = Mi and (5) becomes 
Dy(y)-D-s{y)\ < • {C||/|UA^}^-^ • 

= q.|.c-¥.c'^. 11/11,. AHa-¥). 



Thus we have 



\D"f{y)-D"s{y)\<q^\-C'-^ -C 



(6) 



for < |q;| < / whenever B{y, S) <^ fl for some S > and d is small. 
If d is large. Mi — MqI\5~^ and (5) becomes 



= q,|-Mo-(/!)^-rH 

= q«i-c-||/|U-AH.(z!)¥.ri"i 
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Thus we have 

\Dy{y) - D-s{y)\ < q„| • C ■ (/!)^ ■ r'^l • ||/|U ■ (7) 

for < |q;| < / whenever B[y, 5) C f2 for some 6 > and d is large. 

From (6) and (7) we find that the effect of differentiation appears essen- 
tially only when d is small. 

3 Gaussian Interpolation 

Lemma 3.1 Let h{x) = e"'^''^'^, P > 0, be the Gaussian function in R"'. 
Then given a positive number bo, there are positive constants do, g and G 
which depend on bo and (3 for which the following is true: If f & Ch,m and 
s is the interpolant defined in (1) that interpolates f on a subset X of R"', 
then 

\fix)-s{x)\<A"iGd)HfU (8) 

holds for < d < do and all x in a set Q C i?", where Q is any set in which 
any point x is contained in a subset E of Vt, the set E being a cube or the 
rotation of a cube of side bo satisfying the property that every subcube of E of 
side 2d contains a point of X. Here, \\f\\h denotes the h-norm of f in Ch,m, 
and A" is a constant depending on n. 

Proof. This is just Corollary 1 of [1]. 

Remark: The calculation of A", G and g is quite complicated. The in- 
terested readers are referred to |4j for a detailed discussion. 

Theorem 3.2 Let h, X and s be as in LemmaS.l. For any bo > 0, let do, g 
and G be as in LemmaS.l as well. Suppose Q C i?" and for any x & Q there 
is a set E such that x E E <ZQ where E is a cube or the rotation of a cube of 
side do satisfying the property that every subcube of E of side 2d, < d < do, 
contains a point of X . Let d = d{Q,X) be the fill distance and I > 1 be any 
integer. Then 

(a) for every f G Ch^m o,nd \a\ = I, the distribution D'^f belongs 
to G{R^) and there is a constant Ca such that for all f in Ch^m, 

||D"/-D-s|U<cJ|/|U ; (9) 



6 



(b) for any y E Vt such that the closed hall B(y,5) C Q for some 
S > 0, and any f e Ch,m if ^ <\oi\ <l, then 

\D-f{y) - D'^s{y)\ < q„|Mo"^ • (10) 

where C\a\ is a constant independent of f, Mq — A"{Gd)^\\f\\h 
as in (8), andJii = max{Mi, MqUS'^) with Mi = C"||/||,i where 

Proof, (a) is similar to Theorem2.2 with the only difference that the Gaussian 
function is conditionally positive definite of order m = and Cq, is defined in 
another way because it corresponds to a different radial function h. 

The proof of (b) is also similar to Theorem2.2 with a different definition 
of Mq. We omit it too. (i 

Remcirk: Note that when d is small, Mi — Mi and (10) becomes 

\D-f{y)-D-s{y)\ < q^\{A"{Gd)^\fUy'^.{C'\\fh}^ 

\a\ a /-, \a\ \ 

= C|«|-A"^-^-C"^-||/|U-(G(i)i(i-V). 

Thus we have 

iD^iy) - D-s(|/)| < q,, • A"^-^ • C'^ • ll/iU • (Gd)i(^-^) (11) 

for < |q;| < / whenever B{y, 6) C Q for some 5 > and d is small. 
If d is large, M — Mol\5~'^ and (10) becomes 

\D-f{y)-D^s{y)\ < q^i • Mo'"^ ■ {Mo/!r'}^ 

= q„i-Mo-(/!)^-rH 

= q„i-A"(Gd)i. 11/11,. (/!)¥. rH. 

Thus we have 

\Dy{y) - D"s(y)| < q„| • A" • (/!)¥ • ri"i • • {Gdf^ (12) 

for < |q;| < / whenever B{y, S) C Q for some S > and d is large. 

From (11) and (12) we find that, as in the case of multiquadrics, the effect 
of differentiation appears essentially only when d is small. 
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Finally, we'd like to point out the main goal we are pursuing. In recent 
years mathematicians noticed the importance of the matching of smooth- 
ness. In order to obtain a good approximation, the derivatives should also 
be approximated. This is why we present this paper. 
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